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Abstract 

Recently the 14 moments model of Extended Thermodynamics for dense gases and 
macromolecular fluids has been considered and an exact solution, of the restrictions 
imposed by the entropy principle and that of Galilean relativity, has been obtained 
• through a non relativistic limit. Here we prove uniqueness of the above solution and 

exploit other pertinent conditions such us the convexity of the function h' related to the 
entropy density, the problem of subsystems and the fact that the flux in the conservation 
law of mass must be the moment of order 1 in the conservation law of momentum. Also 
the solution of this last condition is here obtained without using expansions around 
equilibrium. The results present interesting aspects which were not suspected when 
^ ■ only approximated solutions of this problem were known. 

1 Introduction 

The balance equations to describe the 14 moments model of Extended Thermodynamics for dense 
gases and macromolecular fluids are 

d t F + d k F k = o, d t p + d k G kl = o , dtFij + d k G kij = p <tj> , (i) 

dtFm + d k G k m = Pui , dtFau + d k G k uu = Pun , 

where the independent variables are F, F{, Fy, Fm, Fan and are symmetric tensors. P<ij>, Pui, Pun 
are productions and they too are symmetric tensors. The fluxes G k i, G k ij, G k m, G k mi are constitutive 
functions and are symmetric over all indexes, except for k. The restrictions imposed by the entropy 
principle and that of Galilean relativity were firstly studied by Kremer in pQ , [2] , up to second order 
with respect to equilibrium. In [3] we have obtained a non approximated solution through a non 
relativistic limit. To this regard let us remember that the entropy principle for our system (pQ) is 



equivalent to assuming the existence of potentials ti, <p' k and of the Lagrange multipliers A, Aj, A™, 
Xui, X pp u such that 

dti „ dti _ 



0A' 1 d\C 11 dXu 
dti dti 

Fill — TTT > fiill 



dXui ' dXmi ' 
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<9A ' dX{ ' <9AjZ ' 



By comparing ^ 2 with 06 we obtain the following compatibility condition 



b' k _ dti 



dX dX k ' ® 

Moreover, by applying the new methodology [3J, adapted for the present case in [5] and [6], we have 
that the Galilean relativity principle is equivalent to the following two other conditions 

dti dti ( dti dti\ dti 

_ _ Ai + 2A ij — + X jpp ^— 5rs 6 ij + 2— J+4A w?g — (4) 

° " 9A" Ai + 2Xl] dX~ + Xlpp [dXTs ™ * 3 dX~) PPqq dX^i + lk • 
In [3] we have obtained the following solution of eqs. (J3J): 

<f>' k = H V k + H X V? + H 2 V 2 k + H 3 V 3 k , (5) 



^' — 8HqX\ — H1X2 — -H2X3 — -H3X4 



with 



V k = -2X kll (6) 

V\ = —^XkhXhu + ^XppuXk + -XuXku 
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y 2 = —2X k hXhll + -XuXkaXall + ^XkaX a X pp u + 
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~~ ^^ll^kll — XkllXaXail + XkX a iiX a ii + 
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2 12 
+ (t r Kb)^kll ^-XppuXuXk 

V 3 k = 2X pp ii ^2Xl h Xh — trX 2 ab Xk — -XuXk a X a + — Xf t \kj + 

4 17 2 

+ (XkhXh){XallX a ll) — -Xu(X a uX a ll)Xk — —X u XkaX a ll + 

5 25 

— (X a Xaii)XkbXui + {trX 2 ab )XkcXdi + -Xu(X a X a ii)Xku + 
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8 74 4 

+ ^XiiXkhXhii + -^zX n Xku — -Xu(trX ab )Xkii + (XabXaiiXbu)Xk + 
5 o75 5 
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— (X a bX a Xbii)Xkii + -^{trX ab )Xku — 2X k hXhu , 



Xi — Xppu , (7) 
A" 2 = 2A„„A„ JI -yA KJ „A„, 
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X3 — 8X pp u ( — Xf t — -trX 2 b ) + 2X ab X a uX b ii — -XuX a uX a a , 
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X4 = 2X b X a uXbii — tr X b X c uX c ii — -XuXabXaiiXbu + 
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+ 25 ^ll^allX a ii + 8X ppU ( -— A H + -X a (trX ab ) - ^ tr Kb J > 

^5 = --A| + 16A pp zzA-4A a A zz + 2irA 2 6 , 

( 4 1 "\ 

X 6 = 4AA a //A a ;z + 8A pp i/ I — -AAy + -A a A a J + 

8 4 8 4 

+ g A;iA a ^A a — -XutrX ab + yg-^H ~ 4A a feA a AMi + -trX ab , 

X 7 = ^(trX 3 ab )X u - ]^X 2 u trX 2 ab + ^ti + ^X ab X all X bn + 
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+2(trA^)A c A c// - (A a A aZ/ ) 2 - —AXuX a uX a u + 
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+(X a X a ){XbuXbu) — AXl b X a aXb + 

SXppu ^Atr X 2 ab — -X a bX a Xb — — AA^ + — A;/A a A a J + 

12 22 2 

+-^-X[iX a bX a Xbu — —X a X a X a ii . 

X% = ——XuXabXaXbu + 2(trX 2 b )X c dX c Xdu + —XuX 2 b X a Xbu + 
25 5 

148 8 4 

+^Af / A a A aZ/ - - V*rA^)A c A d / + - (tr X ab ) X c X c u - 4X 3 ab X a X bU + 

(8 17 \ 

2X ab X a Xb — (trX cd )X a X a — -XuX a bX a Xb + —X u X a X a \ + 
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+— AA^A a zzA a zz + 16AA pp z/ ( -375^ + ^u( tr ^ab) ~ f rX lb) + 

+^(^)-4^(trA 2 a6 ) + l-^Afz. 

More precisely, our unknown potentials h', <p' k are determined in terms of 4 arbitrary functions 
Hq, Hi, H2, H3 depending on the scalars (0). You can verify that these are solutions of eqs. HD, 
by simple substitution and long calculations. In the next section we will prove uniqueness of this 
solution. In sect. 3 we will impose the further condition ([3|) and, also for this problem, we will find 
an exact solution without using expansions. In sect. 4 we will impose the convexity of the function 
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h! which is important in order that our symmetric system is also hyperbolic. We will find interesting 
results such as the following: Although Xj pp and X pp u are both zero at equilibrium and the first of 
these has an index less than the other, it tends faster to zero when the system tends to equilibrium. 
This fact shows that it is not correct to consider all higher order moments negligible with respect 
to the previous ones. This result confirms the starting point of the new theory called Consistent 
Ordered Extended Thermodynamics (COET) of which we limit ourselves to cite the first paper [7|. 
More precisely, sect. 4 will show that eqs. (|5|) have to be substituted by 
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with Ki arbitrary functions of rji = X\ , rji 
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for i = 2, • • • ,4 and, moreover, of 
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On the other hand, these are compatible with ([5]). Obviously, the form (J8j) can be used only if 

X\ ^ on the initial manifold and until that it remains X\ 7^ 0. 

It is interesting to note that the solution (JSj) , calculated in Xm = 0, becomes 
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On the other hand, if we know (f)' k and h', from the above expression of <ft' k we obtain K\, K2, K3 
because they are coefficients of linearly independent vectors. After that, from the above expression 
of h' we obtain Kq\ also their functional dependence is arbitrary because the above expressions of 
rji - r/8 are the most general possible. In other words, if we know the expressions of (j) lk and b! 
calculated in Xm = 0, we will know them also for Aj;; 7^ ! 

At last, in sect. 5, the problem of subsystems will be considered and, also in this case, we will find 
unexpected results. 



2 Uniqueness of the solution ©-(H). 

In order to prove uniqueness of the solution ([5])- ([7]), let us begin with the case in which the following 
two conditions are satisfied: 

1) The vectors Xm , Xi a X a ii , X 2 a X a u are linearly independent. 

2) The 4-vectors (^T/* 1 ^) ' ( V^j ' ( ' f V^^j 3X6 ^ near ^ independent. 

But, before proving uniqueness of our solution, we need to consider the following representation 
theorem: Every scalar function of our Lagrange multipliers can be expressed as a function of the 
scalars of the set 

Si = {\u , trX 2 , s , trX^, s , X a iiXaii , X ab X a uX b u , X 2 b X a uX b u , X 5 - X$ , X pp u} . 

This theorem can be proved in a way similar to those used for other representation theorems [8], 
PI, [ID], P], as follows: 

It suffices to prove our statement in a particular reference frame and see that in this reference we 
can obtain the Lagrange multipliers from the knowledge of the scalars in Si ; so let us use the frame 
defined by X iU = (X m , 0, 0), A13 = 0, X m > 0, A i2 > 0. 



• If Xm > 0, A12 > 0, we obtain Xm, An, A i2 , A22, A 33 , A 23 respectively from X aU X aU , X ab X all X b u, 
X 2 ab X a iiX b u , X^XaiiXhu, Xu , trX 2 ,^ after that, the 4 th of these can be expressed as function of 
the remaining ones and of trX^. s through the Hamilton-Kayley theorem. 

• If Xm > 0, A12 = 0, with a rotation around the first axis we can select the reference where 
a 2 3 = 0; after that we obtain Xm, An, A22, A 33 respectively from X a uX aU , X ab X a uX b u , Xu , 
trX 2 .^ 

• If Xm = 0, we may select the reference frame where A12 = 0, Ai 3 = 0, A2 3 = 0, and obtain 
An, A22, A 33 from X u , trX 2 rs , trXl s . 

Until now we have obtained Xm and X a b as functions of the elements of Si) obviously, also X pp u is 
a function of them. It remains to obtain A and A&. To this end we note that, from eqs. ((SJ) , ([TJ) it 
follows 




and these are linearly independent for the second hypothesis at the beginning of this section; con- 
sequently the Jacobian determinant, constituted by the derivatives of X^-Xg with respect to A and 



Afc, is non singular. By using the theorem on implicit functions, it follows that we can obtain A and 
Afc in terms of X^-Xg. This completes the proof of our representation theorem. 
So we can now prove our theorem on uniqueness. For the second hypothesis at the beginning of this 
section, we have that it is possible to obtain the scalar functions Hq-H^ such that eqs. © hold. For 
the previous representation theorem, we have that Hi can be expressed as functions of the elements 
in S\. From eq. ©2 we have that also h! satisfies this property, because the coefficients of Hq-H^ 
are proportional to the elements X1-X4 of S\. 

Let us now impose that eqs. ([5]) satisfy eqs. (JH). To this end, let us use the results of [3 J 

dX h dX h ( dX h, c 9X h \ dX h 

— i Aj 1Xij—^\ -|- Ajpp I 7— o rs Oij -\- 2-7-— I + 4App^q 



d\ 1 1 J d\j ' JFP \dx rs rs lJ ' dXijJ ' pmq dx iU 



where h = 1, • • • , 8; r = 0, • • • ,3; Pq = 8, P\ = —1, P2 = — |, -P3 = — |, and there is no summation 
convention over the repeated index r. Consequently, by substituting eqs. ([5]) into eqs. many 
terms give zero contribute and there remain 



dH dH dH 
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with Q 1 = Xu , Q 2 = trX 2 rs , Q 3 = trX^ s . 

Now, for the first hypothesis at the beginning of this section, it follows that the vectors Xm , Xi a X a u, 
Xf a X a u are linearly independent; consequently, the above relation becomes 

This result, for the second hypothesis at the beginning of this section, implies that = 0, that 
is, H r doesn't depend on Qi, Q2, Q3. Consequently it may depend only on Xi-X$, as we desired 
to prove. 

In this way, we have proved uniqueness only if the conditions 1) and 2), at the beginning of this 
section, are satisfied. On the other hand, the set in which these conditions are not satisfied is only 
a sub-manifold of the domain; so our result on uniqueness must hold in any case for continuity 
reasons. This can be clarified better with the following example: If F(x, y) is a continuous function 
such that 



F(x,y) 



5 if y + 
J(x) ify = 0. 

then it follows fix) = F(x, 0) = lim^o F(x, y) = 5 so that F(x, y) = 5 for all values of x, y. 



3 The further condition (US). 



We want now to impose the further condition ([3]); we will see that it can be nicely solved. The solu- 
tion gives Hq, H\, H2, H3, in terms of the arbitrary functions ip = ip(Xi, X2, X3, X±, X$, Yq, I7, Y$), 



t': 



<p = <p{X u X 2 , X 3 , Xt, Z 5 , X 6 , Z 7 , Z 8 ), H* = Ht(X lt X 2 , X 3 , X 4 , Y 6 , Y 7 , Y 8 ) for i going from 1 
to 3, H** = H?*(Xi, X 2 , X 3 , X 4 , Z 5 , Z 7 , Z 8 ) for j = 0,2,3. This solution reads 

*° = i* ( w. + »■ ) 4* (t + »*) + ( w, + + H + < 9 > 



^ , , v f dtp 



T -< 



where it is understood that the right hand sides are calculated in 

Y 6 = X 1 X 6 + lx 2 X 5 , Y 7 = X X X 7 + i X 3 X 5 , Y 8 = XxX 8 + 4 ^4*5 (i ) 
8 12 lb 

Z5 = X 2 X§ + 8X\Xq , Z7 = X2X7 — - X 3 X$ , Z$ = X 2 Xs — - X^Xq . 

In order to prove this result, let us start by noting that from © and (|7|) it follows that V k , V k , 
V k , V k don' t depend on A and, moreover, 
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From dH])9_i2 we have also that the coefficients of Hq, H±, H 2 , H 3 in h! don' t depend on A; 
consequently, eq. (J3J) becomes 



lx 4 ( 2^V k + 2^V k + 2^V k + 2^V k 
2 4 V dX 5 ^ 8X 6 1 ^ 9X7 2 ^ 9X 8 3 
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These equations, for ([TTD 13-16 become 

dX 6 3 dX 7 dX 8 dX 5 3 dX 5 dX 5 ' v ' 

mx dH x 4 dHx dHt mx dH Q 4 dH 2 dH 3 
WXl dX~ 5 ~ 3 X3 dX 7 ~ XA dX, ~ 16Xl dT 6 ~ 3 X3 dT 6 ~ Xi dX G ' 

■t RX dH * 2 X^ 2 X dH2 -^X dHo 2X dHl X 9Hs 

16X1 dx~ 5 ~ 2X2 dx 6 - Xi dx~ 8 ~ 16 1 dx 7 ~ 2 dx 7 ~ XA dx 7 ' 

dH 3 dH 3 4 dH 3 dH dHx 4 dH 2 

16Xl dX~ 5 ~ 2X2 dX 6 ~ 3 X3 dX~ 7 ~ 16Xl dX~s ~ 2 dX~8 ~ 3 X3 dxS ' 

To find the solution of these equations, let us distinguish two cases: 

3.1 The case X x ^ 0. 



1 1 X 4 8H X 1 X 3 dH 2 1 X 4 8H 3 

' 16 X x dX 8 12 Xi dX 6 16 Xi <9X 6 ' 1 ' 
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16 Xi dX 8 8 Xi <9X 7 16 Xi <9X 
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12 Xx dX 7 8 Xi dX 8 12 X x <9X 8 ' 

By substituting these expressions of the derivatives of Ho in (|12|) i . this first equation becomes an 
identity. Let us now change functions and independent variables according to the following relation 

Hi = XxHi(Xx, X 2 , X 3 , X 4 , X5, X1X6 + - X2X5, X1X7 + — X3X5, XiXg + — X4X5) . 

8 12 16 

V V V 

y 6 y 7 y 8 (14) 

for i = 0, • • ■ ,3. With this change, eqs. (fl~3l) become 

1„ ^ 1 „ ~ 1 „ i ' 







0X5 








0X5 








9X 5 





1 



12 16 

— X 3 Ho — — ^ 
12 d 2 16 

— X 3 Ho — — ^ 
12 J 16 



So it will suffice to define ip from 
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X\Hq — —X 2 H\ — —X 3 H 2 — —X4H3 



12 16 dX, 



5 

to obtain, tanks to eqs. (fTlj) . the result but with 93 = 0, H** = 0. On the other hand, from ^ 
we see that the sum of two solutions is still a solution. Consequently, it will suffice now to prove 
that Q is a solution also with ip ^ 0, H** ^= 0, ip = 0, H* = 0; this will be the result of the 
following case. 

3.2 The case X 2 ^ 0. 

Prom eq. (|12P i 34 we obtain 

dffx = dH 2 X 3 dH 1 X 4 ag 2 x 3 ag 2 1 x 4 ag 3 r 5 , 
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dX 8 X 2 dX 5 <9X 6 3 X 2 dX 7 X 2 3X 8 3 X 2 dX 8 ' 



By substituting these in (|12j) -? . this relation becomes an identity. Let us now change functions and 
independent variables according to 

2 1 

Hi = X 2 Hi(Xi, X 2 , X3, X4, X 2 X5 + 8X1X6, Xq, X 2 X-? — - X3X6, X 2 X 8 — - X4X6) . 
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for i = 0, • • • ,3. With this change, eqs. (|15p become 

-8X1^0 + X 2 H\ + -X3^ 2 + -X 4 H 3 

SXiHq + X 2 ^i + -X 3 H 2 + -X 4 H 3 

-8X1^0 + X 2 ^i + -X 3 H 2 + -X 4 H 3 

So it will suffice to define 93 from 





d 


<9X 6 




3if 2 


d 


<9X 6 


dZ 7 


8H 3 


d 


dX 6 


dZ 8 



2 ~ 1 dip 
-8X\Hq + X 2 Hi + —X 3 H 2 + — X 4 H 3 
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to obtain, thanks to eqs. (fT6|) . the eqs. Q, but with ^ = 0, iT* = 0, as afore said. 



4 The convexity of h! . 

In order that our system ([1]) be hyperbolic, we have now to impose that the hessian matrix Qx A d\ B 
is positive defined, with Aa the generic component of the Lagrange multipliers. In other words, 
the quadratic form Q = Qx A d\ B S^aSXb has to be positive definite. Let us exploit this with the 
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potentials ©i; in these expressions, except for replacing Xi with X; L /(X\) for i = 5, • • • , 8, the 
remaining polynomials in Xj/ (X\) for j = 2, • • • ,4 have been chosen in order to eliminate from Xi 
for i = 5, • • ■ ,8 the terms depending only on X au . 

Well, we want now to evaluate this quadratic form Q in the state, which will be called C, where 
Aj = 0, Ajj = \Xu5ij, Xui = 0; so there remain, as independent variables A, Xu, X pp u. This is 
an intermediate state with respect to equilibrium, where we have also X pp u = 0. To this end we 
need the expressions of our variables up to second order with respect to the state C. After some 
calculations, we find 
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with A <a b> = X a b — ^XuSij. 

After that, from h' = h'(rn), we find that the expression of h' up to second order with respect to 
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the state C is 
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where the apex * denotes a quantity calculated in the state C, so that we have also 

* v \ 16 x * 32 \2 * 64 \3 

rjl = Xi = A pp u , r) 2 = — —An , r) 3 = — A u , r? 4 = - 12 5 " ' 

* ic\ * 32 \ \ * 64 \\2 * 8-16 3 
7] 5 = 16A , t/6 = — g-AA« , rj 7 = —XX n , % = " 27. 125 

Taking into account these intermediate results and the additivity of Q = q\ a q\ b SXa$Xb and by 
calculating g^g Xg in the confront state C, we find 

Q = Qi + Q2 + Qs , with 
Qi = aii(5A) 2 + 2ai2<5A<5A;/ + 2ai3(5A<5A pp ii + a 2 2(5An) 2 + 2a2^Xii5Xp P u + a^(5\ P pu) 2 , 

Q 2 = 6 n ( p!L ). ( 6 -^L) + 26 12 ( ^) *A a + b 22 SXJX a , 

Q3 = c <5A< rs ><5A< rs > , 
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d 2 h'(yt) 
dXudXppu 
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Consequently, the required convexity holds if 
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We note that these conditions are continuous in X pp u, so that we may impose them only calculated 
in Xppu = 0; in this way we will obtain the requested convexity not only in a neighborhood of the 
state C, but also in a neighborhood of equilibrium. 

We have performed the same passages also by starting from eqs. ©i, instead of ©i; in this way 
we have found that Q isn't positive defined. We conclude that only eq. ([8]) is the correct expression 
to use. 

We note that also the results of the previous section can be written taking into account the expression 
©i. In particular, we can use the expressions at the end of section 1 to find X\-Xg as functions of 
771-778. After that, by using also eqs. (fl"0]) . we can obtain Y^-Yg, with Y§ = X5, that is 



Y 5 
Y 6 
Y 7 
Y* 



VlV5 - fjT/3 + ^(%) 2 



(17) 



m 
m 



11 1 / x2 



1 1 , x2 1 



24 



1 

16 



7?1778 + ^771774775 + —^774 (77 2 ) 2 



1 

2io' 



32 



773774 



Prom Ki = rjiHi and by defining # = rjii/j, K* = ^H* for 7 = 1,2, 3, we can rewrite eqs. ([9]). We will 
limit ourselves to the case X\ ^ 0, so that we have <p = 0, JET?* = 0. The result is that the solution 



1 o 



gives Kq, Ki, K2, K3, in terms of the arbitrary functions 1? = #(771, 772, 773, 7/4, Y5, Ye, 5^7 5 ^s)> 
if* = K*(r]i, 7/2, 7/3, 7/4, lg, I7, Yg) for i going from 1 to 3. This solution reads 



1 (d-d A 1 /St? A 1 / 00 A 90 



dY f 



K * = m[— + K* 3l , 



where it is understood that the right hand sides are calculated in (fTTj) . 



5 The subsystems 

Other interesting particulars of our solution become manifest when we search the subsystems of ([I]). 
As example, eqs. (j3|) calculated in Xm = 0, X pp u = become the conditions we would have by 
starting only with (JTJ) i_3 . But eq. ©2 in Xm = 0, X pp u = gives h' = which cannot be 
accepted for the required convexity. This problem isn't avoided neither by using eqs. ([8]) because 
the consequent solutions don't satisfy the conditions ([!]) calculated for the subsystem. To verify 
that this is the case, it suffices to note that (jl|)i with 775 instead of h! is satisfied, but if we replace 
775 with its value in Xm = 0, X pp u = 0, that is 16A, we see that this satisfies no more eq. ^3J) 1 
calculated in Aj^ = 0, X pp u = 0! The reason is evident from the fact that 775 satisfies eq. (@J)i; but, 
if we calculate this equation in Xm = 0, we find 

n dri5 \ -lo\ dr l5, A , ( d7? 5 \ 



= — A, + 2A ii — -16A,, 

whose value in X pp u = isn't a solution of eq. (jl])i calculated in Xm = 0, X pp u = 0. 

An idea may be that to redo the passages of section 5 of paper ([3]) but starting from the beginning 

with Xm = 0, X pp u = 0, that is, with 



\P — _L 

m 



X 13 = — 
m 



'f Ajj 



0; 



Xij + hXuSi 



-§Au 

0; 







Ml 









1 

c z 



-X 



0, 



but in this case, among the scalars there is the one coming from XpX^ +^Q\rn^c 2 , that is A a A a — |AA// 
and this, substituted to h' in (j3J)i calculated in Xm = 0, X pp u = 0, doesn't satisfy it. The same thing 
can be said if we start from eq. (89) instead of (80), both of (|_3j). While, if we start from eq. 



1 o 



instead of (80) (both of (|3j)), we obtain quickly h' = 0, 4>' k = 0. In other words, the subsystem 
with 10 moments cannot be obtained in any way as a non relativistic limit. 
What about the subsystem with 5 moments? 

If in eq.(9) of ([3]) we substitute Ay = |An<% , Xm = 0, X pp u = 0, they become the entropy principle 
for the system constituted only by ([11)1,2 and by the trace of ([1])3, with Lagrange multipliers A, Aj, 
^Xu respectively. 
The equations then become 

dh' 2, dti d(f>' k x 2, d(j)' k , /c 
oX 3 oXi oX 3 oXi 

With arguments like those above described, the solution of this equation cannot be found from KSJ), 
nor from (jHJ) calculated in the above values of Ay, Xm, X pp u. 

Instead of this, the idea of redoing the passages of section 5 of ([3]), but starting from the beginning 
with Ay = ^XuSij, Xm = 0, X pp u = 0, is successful . In fact, starting from (80) or from (89) (both 
of (PI)), we find 



h' = -h u H , <j,' k = H Q X k 



3 

where Hq is a function of Xu, X a X a — |AA^. These functions satisfy effectively eqs. (fT8|) . More than 
that, we have that they satisfy automatically also eq. ([3])! 

Obviously, we cannot obtain this result by starting from the beginning from (88) of (|3j) because in 
this case we would obtain quickly h! = 0, (f)' k = 0. On the other hand, if we start from (88) of ([3]) 
we obtain A' 3 = and this isn't adequate to describe the relativistic model; the less for its limit! 
In other words, we have found that the following diagram isn't commutative 

(relativistic model — > — > — > (relativistic subsystem with 
with 14 moments) with 5 moments) 

J, non relativistic limit J, non relativistic limit 

(classical model (classical subsystem (classical model 

with 14 moments) — ► with 5 moments) ^ with 5 moments) 

and this is quite different from the results obtained with expansions around equilibrium, that is, 



(classical model (classical subsystem (classical model 

with 14 moments) — > with 5 moments) = with 5 moments) 

even if this has been until now proved only for the less restrictive case of ideal gases [12 
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